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Abstract 

We calculate the reduced density matrix for a system of coupled harmonic 
oscillators in a bosonic heat bath using the Born-Markov approximation and 
show that the expectation values of position and momentum evolve like classi- 
cal quantities. We consider the cases where every oscillator is coupled to its own 
reservoir and where all oscillators are coupled to the same heat bath, pointing 
out the differences between these two models. In particular, for systems cou- 
pled to a common reservoir, we present an uncommon dissipation mechanism, 
which disappears, when the dynamics of the system are gouverned by the low 
frequency modes of the environment (Caldeira-Leggett limit). We also address 
the interesting phenomenon of reservoir-induced interactions. 

Keywords: decoherence, Quantum Brownian Motion, Born-Markov 
approximation, interacting quantum systems 



1. Introduction 

In the early 1980s decoherence theory made a first attempt to explain the 
emergence of classicality in quantum mechanics [l|. While the quantum-to- 
classical transition is today well understood [2|, y, |J, |5j , decoherence is nowadays 
the major obstacle in the experimental implementation of quantum information 
processing [f|. It is therefore a challenging task to find mechanisms that sup- 
press decoherence. Recently, dissipation and decoherence free subspaces have 
been found in oscillator networks yj . The underlying mechanism is that in mul- 
tipartite systems interactions between the subsystems can lead to the robustness 
of certain states, while other states remain affected by decoherence. A second 
reason why one is interested in the decoherence of interacting quantum systems 
is quantum synchronization. In classical mechanics, two coupled dissipative sys- 
tems can evolve into a synchronized state [8] and will stay there for all times. 
If one wants to find a quantum analogue of synchronization, one is confronted 
with the problem that quantum mechanics is invariant by time reversal, which 
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provides no reason why two initially synchronized quantum systems could not 
desynchronize again, simply by reversing the arrow of time. It is suggested that 
decoherence, i. e. the interaction of a system with its environment, as it breaks 
time reversal symmetry, is an essential ingredient in order to study quantum 
synchronization [t| ■ 

As the interaction of a system with its environment can be very compli- 
cated, it is clear that analytical calculations are limited to a very small class 
of models. One of the simplest models is a harmonic oscillator linearly coupled 
to a bosonic heat bath and has been studied ingreat detail in literature under 



the name of Quantum Brownian Motion ll|, Il2l Il3l . Il4l 1 1 51 ] . Although some 
work has already been done on the decoherence of interacting quantum systems 
0, Gjl EH ! the problem of several interacting harmonic oscillators is not yet 
fully understood. As an important contribution, Cacheffo et al. [l9| solved the 
problem in the so-called Caldeira-Leggett limit. While these authors use path 
integral methods to derive the reduced density matrix [20], in the present work 
we will apply the Born-Markov approximation, which not only reproduces the 
results of [19j but also gives more general results and is less complicated from 
a mathematical point of view. In the present paper, we will derive and solve 
the Born-Markov equation of N linearly coupled harmonic oscillators, each one 
interacting with its own bosonic heat bath, and point out the formal analogy 
with the problem of a single oscillator in a heat bath (section^. After that, we 
will compare the result for distinct reservoirs with that obtained for the coupling 
to a common reservoir (section [3]) . In detail, we will consider a general spectral 
density and analyze the behaviour of the expectation values of position for the 
case of two interacting oscillators (section . 



2. Distinct Reservoirs 

In the present section we consider the system Hamiltonian 

Hs = J2 (oS- + \ M ^l X l) +\H 9apX a X , (1) 

a=l ^ Q ' a ^/3 

where X a and P a are the position and momentum operators of N harmonic 
oscillators with frequencies Q a and masses M Q , linearly coupled by coupling 
constants g a p — gp a . Every oscillator a is coupled to its own bosonic heat bath 




a) a 2 a 2 
{a) 2 3 3 3 



(2) 



(where u>^ are the frequencies and the masses of N^g non-interacting 

harmonic oscillators) via an interaction Hamiltonian 



N 



(a) 



H\ a) =X a ®Y. ^X a ®E a (3) 
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with coupling constants c " . The total Hamiltonian is therefore: 



(a) 

H = H S +H £ + H I (4) 



with 

N N 



2.1. Derivation of Born-Markov Equation 

In order to determine the temporal evolution of the system described by ^ , 
it is useful to have a master equation for the reduced density matrix ps of the 
system. In the case of Quantum Brownian Motion an exact master equation 



has been derived by Hu, Paz and Zhang [14J. A much simpler form can be 
obtained by applying the Born-Markov approximation, which means that the 
system-environment interaction is weak enough so that the full density matrix p 
of the system and the environment approximately factorizes for all times (Born 
approximation) : 

p(t)*ps(t)®p e (Q) (6) 

(where ps(0) is the initial state of the environment) and that the correlations 
(EcP(t)Ejp (t')) pe of the environment decay more rapidly than p$ changes 
(Markov approximation) . Here, {■) P£ denotes the expectation value with respect 
to pf(0), while the index (J) denotes the temporal evolution in the interaction 
picture, i. e. 

E^\r)=e^ Hs+H ^ T E a e-^ Hs+H ^ T . (7) 

Using the Born-Markov approximation, one can derive the following master 
equation [2j: 

J POO 

- Ps {t) = -i \H Sl ps(t)\ -J2 dT (C af >(T)[X a ,XjP{-T)p S (t)] 

M <*,p J ° (8) 

+C af} (-T)[ PS (t)Xj i I \-T),X a ]) , 

where 

C aP {T):=(E<P{T)EjP(Q)) (9) 

N 1 PS 

is the correlation function between the baths at time difference t. 

Expressing Xj and in terms of annihilation and creation operators 

{xf ] = {bf ] + bf')/fi^Jf, p\ a) = iy/m^ a) /2(b^ - bf)) and 
using the Baker-Campbell-Hausdorff formula 

oo 1 

e x Ye- x = J2-[X,Y] n (10) 

n=0 
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(with the nested commutator [X, Y] n := [X, [X, Y] n _i], [X, Y]o := Y), one finds 
for the correlation function 



(11) 



Inserting the equilibrium density matrix 



g -n £ /Jo 

(where T Q is the temperature of the a-th heat bath, k B = 1) and evaluating the 
trace in the occupation number representation, one immediately sees that 

bfb^) =(bfb[^) = (13) 

I PS \ / P£ 

' pe e " 3 /r Q _ i \ I P e y g^- / T = _ i / 

(14) 

(j = 1, . . . , ivj£g and I = 1, . . . , Nffg). The correlation function can then be 
written as 

C a f}{T) = S al3 (v a (T) - ir) a (T)) (15) 

with the noise kernel 

n hb (a) 2 / (<*)\ 

Mr) - g coth ^ j coe(a,<">r) (16) 

and the dissipation kernel 

<°2 c («) 2 

Calculating X„\t) is a bit more complicated. Again it is easier to work 
with annihilation and creation operators a a , a^ a (a = 1, . . . , N). Observing that 

N 

[H s ,a a ] = -Sl a a a - ^ 9ap( a P + a p) ( i8 ) 

(with 5 Q/ 3 = g a p / y/4M a Sl a Mp£lp), one is led to 

[H s ,a a ] n = h^ a, (19) 
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where a = (a\, a%, . . . , ajy, (in) t and h„ is a vector valued function depending 
on fl a and g a p (a, = 1, . . . , N). By iterating (fT9|) one finds that 

[H s ,a a ] n = G n h<, a) -a, 



(20) 



where the matrix G is defined by 

/Wi r 12 

G := Fl2 



\TlN 

n a o 
o -o Q 



T 1N \ 



rjv-i,jv 



■ Q/3 



(21) 



(22) 



and (ho"'') — ^i,2a- It follows that the evolution of X a is given by 



N 



0=1 



X P\I || II ' (cOs(Gr) 2/ 3-l,2a + COs(Gr) 2 /3,2 Q ) 



ft 



(sin(Gr) 2/ 3-i,2 a - sin(Gr) 2 ^ i2Q ) 



(23) 



where cos(Gr) Q ,. / 3, sin(G7") a)( g are matrix elements of the matrix valued cos/sin 
function. 

The Born-Markov master equation (|5J| can then be written as 
d N 

-p s (t) = -i[H s ,ps(t)}- (Dap[X a) [Xp,ps(t)]] + i7 a p[X a ,{Pp,ps(t)}] 

a,0=l 

+f a p[X a) [Pp,Ps{t)}} + l -M^[X a ,{Xp,p s {t)}\ 



with the following coefficients: 



O z — 

la0 ■ = 

D a p := 



wj dT ^ T) 

dr ?7 q (t) 



y/MJlaMpfip 



(24) 



(cos(Gr) 2/3 _i j2a + cos(Gt) 2/3j2q ) (25) 



(sin(Gr) 2/ 3_i i2ce - sin(Gr) 2 ^ ;2a ) (26) 



<lr (>)\/ -jtf-Q - (cos(Gt)2^i,2„ + cos(Gr) 2/3;2a ) 



(27) 



/•°° 1 
fa0-=~ &Tv a (T)— ====== (sin(Gr) 2(3 _i !2Q - sin(GT) 2/3 , 2Q ) (28) 
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2.2. Interpretation of Born-Markov Equation 

Let us now compare the Born-Markov equation (|24|) for coupled oscillators 
with the one for a single oscillator, which states [2|: 



d^ (t) = 



^ + l M (n 2 + n 2 )x 2 , Ps (t) 
(D[X, [X, p s {t)]] + n[X, {P, Ps (t)}} + f[X, [P, Ps (t)}]) , 



where M and £1 are the mass and frequency of the oscillator and 

i 

dr i](t) cos(fir) 



2 

M 
1 



7 == 



dr 7/(r) sin(57r) 



D := 



/:= 



dr v(t) cos(fir) 



MQ 



dr v(t) sin(f2r) 



(29) 

(30) 
(31) 
(32) 
(33) 



It is clear that the coupling of the oscillators between each other leads at the 
level of the master equation to double commutators containing operators that 
belong to different oscillators. Furthermore, the coefficients (f2"5]) -([2"5 f for the 
coupled oscillators are no more Fourier transforms of v a {r), r/air) with respect 
to the bare frequencies fl a but with respect to the eigenfrequencies of the matrix 
G. In the limit of vanishing coupling g a p — ¥ the matrix G becomes diagonal, 
which leads to = D a p = j a p = / Q( g = for a 7^ j3. As expected, one 
obtains N independent decohering oscillators with coefficients (f3"0)) - (|c>3"l) . 

In reality one is often confronted with a continuous distribution of frequencies 
and coupling constants of the environment. In this case one can introduce 
spectral densities J a (ui) (which together with the temperatures T a describe all 
aspects of the reservoirs which are relevant for decoherence) , so that (fT6|) and 
(fTTj) write 

v a {r) := J duj J a (uj) coth (^r^j cos(wr) (34) 

/•OO 

^o,(r) := / dui J a {uj) sin(wr). (35) 
Jo 

At this point, a comparison with the master equation obtained by Cacheffo et al. 
in 19] is interesting. The authors of [lj} consider the Hamiltonian (j4]) for a — 2 
oscillators but with additional position-momentum and momentum-momentum 
couplings and apply the Caldeira-Leggett approximation. If in (f24| one uses 
ohmic spectral densities (see Q for technical details and [ll| for the physical 
meaning of ohmic dissipation) 



(36) 



6 



with cut-off frequencies A a and assumes that the typical frequency scale f2 = 
max(fii, SI2) of the system is lower than the frequency scale A = min(Ai, A2) 
of the environment and that the temperatures T a are high compared to the 
frequencies of the environment (A Q <T a ), one obtains the master equation 

±ps(t)=-i[H s ,ps(t)} 



2 

£ 



2M al { a) T a [X ai [X a ,p s (t)]\ + i^ a) [X a , {P a ,p s (t)}] (37) 



iM a ^ a) A a [X a ,{X a , Ps (t)}} 



derived in |l9j for two linearly coupled oscillators, provided that the position- 
momentum and momentum-momentum couplings introduced therein are ne- 
glected. As the authors of [3] remark, the Caldeira-Leggett limit with only 
position-position coupling is trivial: it generates no new non-unitary contri- 
butions. It is interesting that the master equation (38) in [Til ] with (bilinear) 
position-momentum and momentum-momentum couplings is very similar to our 
master equation (|24|) (of course the terms in f a p are absent, as they vanish in the 
high temperature limit). Note that there is no problem to introduce position- 
momentum and momentum- momentum couplings in our model as well, this 
would merely change the matrix G. As long as the Born-Markov approxima- 
tion holds, we have therefore extended the master equation from (l9j to arbitrary 
temperatures and spectral densities. 

To see the physical meaning of the coefficients (l25l) - ([28|) . we set up differen- 
tial equations for the expectation values of position and momentum. This can be 
done using d/dt(0)t = Tr(dps/diO) for time- independent observables O and 
inserting (I%4"1) . In this way one obtains the following Ehrenfest-like theorem: 

d /(X)A _ f(X)S 
dt l(P)J " Adiss l(P> 



(38) 



where the position and momentum operators are contained in the vectors X :— 



{Xx, 
by 



with 



, Xn) an d P := (Pij • • • 7 Pn) and a dissipation matrix Adi SS is defined 



Adiss ■ — 





-g' 



M 



-27 



M :=diag(Afi,...,M 



N 



521 



V 9' m 



7 == 



512 



711 



\7ivi 



9n,n-i 



71 AT 
7-/V7V 

9in \ 

9n-i,n 



(39) 



(40) 



(41) 
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Here we defined the renormalized coupling constants g' aB :— g a p + MpQ^g 
(note that g' aB ^ g' Ba in general) and the renormalized frequencies := fi^ + 
fl^a ■ The coefficients ti 2 aB therefore merely renormalize coupling constants and 
oscillator frequencies. From t|38[) one can see that "/ aa describes the dissipation 
of the momentum of the oscillator a and that (a 7^ f3) is responsible for the 
dissipation of the oscillator a caused by the coupling to the oscillator (3. The 
formal solution of (15"8l is: 



! X )t \ _ A diss t ( (x) t=0 

;p>J~ e l(P) t =o 



(42) 



The differential equations for the second moments of X a and P a - in contrast 
to (|38[) - contain also the temperature dependent coefficients D a g and f a B- 

±(X a X 8 ) t = ^({X a ,P,}h + ^- a ({XM) t (43) 
d N 

jA p a Pp)t =D a p + D Ba - 2j2(l a s(P a Ps)t + lps{PpP&)t) 

5=1 

- ^M B ^{{X B ,P a }) t - ±M a n'*({X a , P B }) t (44) 

N N 

E 9 , l)S {{x s ,p«})t+ E ^({^.-p/j})* 

I .... ., . ^ 2 



-({X Q ,P^}) t = - 2/^ ^^^P^*}}, + j^(P a Pp)t 

S—l 



(45) 



AT 

-2M 3 nf{X a X 3 ) t -2 E fl^W^t 

As in the case of Quantum Brownian Motion Q , the coefficients D aa describe 
diffusion in momentum and decoherence in position basis of the oscillator a. 
Therefore, the "off-diagonal" coefficients D a s account for the decoherence of the 
oscillator a caused by the coupling to the oscillator j3. A similar interpretation 
can be attributed to the anomalous diffusion coefficients f a g. 

2.3. Solution of Born-Markov Equation 

In contrast to the exact master equation, the Born-Markov equation (|24[) 
can be solved easily. To do this, it is useful to change to a more suitable basis 
via the transformation [3, HI 

p(k, A, t) := Tr (exp (*(k ■ X + A ■ P)) Ps (t)) = Tr (e ikA /V k ' x e iA - p p s (t)) , 

(46) 

where k := (fci, . . . , fejv), A := (Ai, . . . , Ajy). Physically, the coordinate A 
measures the distance from the diagonal x = x' in position basis, while k can 
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be associated with the pos ition on the diagonal. Note also the formal similarity 
to the Weyl transform [21] . By applying (|46|) on (|24)l one obtains the following 
first order linear partial differential equation: 



d_ 

di 



p(k,A,t) 



N 

E 



a=l 



a <9A a <% Q <9A a 



p(k,A,t) 
(47) 



Equation (|4"T)) can be solved using the method of characteristics. On a char- 
acteristic line (k(r), A(r), t(r)) the original equation reduces to an ordinary 
differential equation: 



-p(k(r),A(r),t(r)) = (k(r), A(r))B doc (k(r), A(r)) T , 



where the decoherence matrix Bd cc can be defined by 

Bdcc '■ — 
( fll ■ ■ ■ flN 





-f D 



/-Dn 



D 



V/- 



Nl 



NNj 



D 



IN 



(48) 



(49) 



(50) 



\£>jvi 



D 



NN / 



By comparing 07} with ^ = ^ • V k p + ^ • V A p + one obtains: 



dr dr \ A ( T ! 

Solving (|5l"j) and inserting in (|48|) yields: 

P(r) 



-^-diss 



k(r) 
A(r) 



(51) 



p(k , A , 0) exp f(k , A ) jf e Ad - r 'Bdcce A ^ T 'dr' (k , A ) T ^) (52) 



with k(r = 0) =: ko, A(r = 0) =: Ao, i(r = 0) = 0. After retransformation 
into the time domain this gives the compact form 

P (k, A, t) = Po ((k, A) e A ^<) exp (-(k, A) C(f)(k, A) T ) (53) 

with initial condition p(k, A, 0) =: po(k, A) and decoherence kernel 



ft 

C(t) := e Adi3S * / e - Adi « T B doc e - A d, s ^ dr e A * 
Jo 



(54) 
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The (k, A)-representation is related to the position representation via 

(x|p s (t)|xO - J d N k e- lk ^p(k, x - x', t) (55) 

(with x := (xi, . . . , xm) t : x' := (xj, . . . , x' N ) T ), as can be verified by inserting 
in (1461) and evaluating the trace in the position basis. 

It is interesting from a mathematical point of view that the solution ([55} 
factorizes into a part depending on the initial condition and a Gaussian part 
describing decoherence and therefore independent of the initial condition (a sim- 
ilar result was found for Quantum Brownian Motion in [22| ) . It is also surprising 
that the same matrix Adi ss that describes the evolution of the expectation values 
(X) and (P) enters in the part of the solution depending on the initial condition. 
If one compares (|53|) with the solution of the single oscillator problem, one finds 
that the two solutions are formally equivalent. The only difference is that in 
the case of the single oscillator, the dissipation matrix Adi ss and decoherence 
matrix Bd ec are defined by: 

.(iv=i) ._ / 1/AA (JV =i) ._ ( 

A diss — \-MQ, a -2 7 J doc — \-f D 

As expected, in the case of the coupled oscillators mixed terms k a kp, k a 
and A a Ap (a ^ j3) appear in the Gaussian. In order to explicitly calculate the 
time-dependent matrix C(t), one has to diagonalize Adi SS - Interestingly, this 
can be done analytically only in the case of N = 2 interacting oscillators (for 
higher N the order of the characteristic polynomial of Adi ss will be greater than 
4 and therefore cannot be solved analytically, if no additional symmetries are 
assumed). 




3. Common Reservoir 

Now, we consider the case where the system described by the Hamiltonian 
H$ is coupled to a common reservoir 

Nhb ( 2 . \ 

via an interaction Hamiltonian 

N N HB N 

= x ° ® E c f x i ^Y. x «® E « c) - ( 58 ) 

a— 1 J — 1 a— 1 

The time evolution of the operator X a in the interaction picture is therefore 
still given by (|23p. but the correlation function is now different. Using the 
annihilation and creation operators bj, &t (j = 1, . . . , Nhb) of the heat bath, 
one obtains 

Nhb J a )J^ . . w . . 

Cil (t) = V . J ' ( ( 6je^ T + 6,-c-* J i T ) (6/ + &, f ) ) . (59) 
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Consequently 

C${T)= V ${T)-iT$(T) (60) 

with noise kernel 

N HBc (a) c (p) oo 
p e$( T ) : = H 2m L C ° th V2t) cos ^' r ) = J duj coth (27^) cos(wr) 

(61) 



and dissipation kernel 

WffB „(<*) J/3) 



»?ifl( T ) : =S-5 — — sin(c^T)= / dw J Q(3 (o;)sin(wr) (62) 
3 - =1 * m jU Jo 

(where T is the temperature and J a p(u) the spectral density of the heat bath). 
Comparing (|60|) with (flST) . one sees that in the case of a common reservoir 
the correlation function is no longer diagonal in the indices a and j3. In the 
following, we will only consider the case of equal coupling of the oscillators to 
the heat bath, i. e. Cj = cf = Cj and consequently J a p(w) = J{w) (for all 
a, (3 = 1,,.,,N). Using (|60|) . one is therefore led to a different Born-Markov 
equation: 

d N 

-^Ps = -i[H s ,ps\- X! ( D p[ x a>[ x l3,Ps]]+i7p[X a ,{Pp,ps}] 

a, 0=1 (63) 

+fp[X a , [Pp,ps]] + l -M^}[X a ,{X^p s }] 

where the coefficients are sums of the coefficients (j2"5)) - (j2"5)) for distinct reservoirs, 
i. e. 

N 



Y := J2 Y *0 (64) 



with Y = O 2 , D, 7, /, taking T a = T, J a (u>) = J(w) Va = 1, . . . , N . Interest- 
ingly, the structure of the Born-Markov equation for a common reservoir is the 
same as for distinct reservoirs, only the coefficients, which contain information 
about the influence of the heat bath on the system, are different. Therefore, 
the density matrix for distinct reservoirs is formally equivalent to the one for a 
common reservoir. Physically, (|64[) means that the oscillators are not only cou- 
pled by an explicit interaction but also indirectly via the interaction with the 
heat bath (figure [T]). This feature of the common reservoir becomes particularly 
striking, if one considers the limit of vanishing coupling g a p — > 0. In this case, 
one obtains the master equation 



9 Ps = -i[H s ,p s ]- J2 (D^ 0) [X a ,[Xp,p s }}+i7 l r 0) [ x a ,{P^Ps}] 

a, 0—1 

+f^= 0) [X a , [P ,ps]} + 1 -M P (f^ =0) ) 2 [X a ,{Xp,p s }} 



dt 

a, 0—1 



(65) 
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Despite the vanishing coupling constants g a p, the oscillators are mutually cou- 
pled by a reservoir-induced interaction, which is due to the fact that all oscilla- 
tors are coupled to the same reservoir. In the next section, we will investigate 
this coupling mechanism in more detail. 



N 




heat 


y 


bath 




Figure 1: Illustration of the influence of a common reservoir: The coupling to the same 
reservoir introduces reservoir-induced interactions, therefore e.g. the dissipation coefficient 71 
of oscillator 1 is the sum of the dissipation coefficients j a i for distinct reservoirs. 



4. Case of Two Interacting Oscillators 

4-1- Dependence of the Coefficients on the Interaction 

In this section, we will consider a specific physical situation described by the 



commonly used spectral density |12l Il4 123| 



J(u) = -Mfow'A^'e-"'*, (66) 

where A is the cut-off frequency and 70 > the coupling strength between the 
system and the reservoir. The exponent s > describes the type of dissipation: 
< s < 1 sub-ohmic, s — 1 ohmic and s > 1 super-ohmic dissipation. By 
varying the exponent s, one is a able to see the passage from local (s = 1) 
to non-local dissipation (s 7^ 1) and from white (s = 1, high temperature) to 



coloured noise |12|, |14| . 



If one wants to understand the physics of the system, the first step is to 
evaluate the coefficients appearing in the master equation. In the following, we 
will concentrate on the two-oscillator system with equal masses M\ = Mi 
\ 1 and coupling constant 512 =521 =9 and will assume the same spectral 
density (|66[) and temperature T for both reservoirs. Furthermore, we will focus 
on the coefficients 7 a/ g and D a p, which are particularly easy to calculate, as 
they are double Fourier transforms of J(w) and J(u) coth(w/(2T)), respectively. 
Moreover, 7 aj g and D a p are expected to control the dynamics of decoherence, 
as it is the case with 7 and D for Quantum Brownian Motion [4|]. For distinct 
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reservoirs one obtains: 

7o 

7n = 



2Jk 2 +r 2 



k 2 + r 2 _ + r_ j (r + + yj ' k 2 + r 2 _ 



[ yj 1 n 2 +r 2 _ -f-j ( r +- \J n 2 + r 2 _ 



(s-l)/2 



(s-l)/2 



(67) 



712 = 721 



70 K 



2J K 2 +r 2 



\J k 2 + r 



-i)/2 , == 



. \ («-l)/2 , = 

- ^k 2 + r 2 ) e -V r +-V^+' J 



(68) 



n 



7o MA 
2,/k 2 +r 2 



- r- ; -f /■_ ](/■_ + y'ft 2 _|_ v 



x coth 



s/2 



.1 i ~2 I „ \ / „ I 



2T 



V 



(69) 



/ 



+ (^^J k 2 + r 2 _ - r^j - \J k 2 + r 2 ^j 



s/2 

+ r±) e 



:- j_ ,-2 | V /, '+ _ V" K2+r - 



x coth 



■+ — y k 2 + r 2 _ ^ 
2T 



#12 = #21 



7o«MA 
2^K 2 +r 2 _ 



s/2 

r + + \l k 2 + r 2 _\ e 



+rl 



coth 



- (r+ - \J k 2 + r 2 _^j 



s/2 



coth 









2T 




J 




- ^Jk 2 + r 2 ^ 


2T 



(70) 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

g/(MA 2 ) g/(MA 2 ) 

Figure 2: Dependence of 711 (left) and 712 (right), normalized by 70, on the coupling strength g 
(here: r+ = 1, r_ = 0). The dotted curves correspond to s = 0.5 (sub-ohmic), the continuous 
ones to s = 1 (ohmic) and the dashed ones to s = 2 (super-ohmic). For subohmic dissipation 
7n (712) diverges strictly monotonically to +00 (—00), when g — > g^. In the case of ohmic 
dissipation, one finds the same monotic behaviour, but no divergence occurs. For super-ohmic 
dissipation there are different types of curves, ranging from strictly monotonically increasing 
and decreasing curves to curves having a local minimum or maximum (here only the case of 
a strictly monotonically decreasing/increasing curve is depicted). 



where r± := (ft? ± Vt 2 )/{2k 2 ), k := g/(MA 2 ), f := T/A. The coefficients 
722 and D22 can be obtained from 711 and D\\ by making the substitution 
r_ — > —T-. Before starting to analyze 7 Q) g and D a p, it should be said that all 
coefficients in the master equation are only defined for g 2 < M 2 tt 2 £l 2 =: g 2 ^. 
This is an anomality due to the specific choice of the coupling between the os- 
cillators, which occurs when one diagonalizes H$. 

If one takes a closer look at the dissipation coefficients 7 Q/ 3 (figure [5]), one 
sees that for sub-ohmic dissipation they are strictly monotonic with g and di- 
verge in the limit g —¥ goo (to +00 for a — f3 and to —00 for a 7^ /3). For 
ohmic dissipation 7 Q ^ is also strictly monotic with g (increasing for a = /3 and 
decreasing for a ^ (3) but does not diverge. It is also interesting that for a ^ (3 
the dissipation coefficients are always negative or zero in the case of ohmic and 
sub-ohmic dissipation. On the contrary, for a = f3 the coefficients are always 
positive, independently of s. In the super-ohmic case, the form of 7 Q( 3 as a 
function of g explicitly depends on r± and T. In particular, 7^/3(3) can present 
a minimum or maximum for values of s near 1 . 
Concerning the decoherence coefficients D a p 1 one finds that for high tempera- 
tures T 3> max(J7i, f^) the relation 

D aP » 2MT laP (71) 



holds. This relation is well known from the Caldeira-Leggett equation [llj . but 
in contrast to this, (fTTj) not only holds for ohmic but also for sub- and super- 
ohmic spectral densities. Moreover, the above result is completely independent 
of the cut-off frequency A of the environment, which is not the case for the 
Caldeira Leggett approximation (see section 12.21) . Consequently, for high tem- 
peratures, D a p depends in exactly the same manner on the interaction strength 
g as 7 0j g, the temperature simply rescales the D Q /3-axis (see figure El only flu 
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Figure 3: Dependence of Dn, normalized by ■ygMA, on the coupling strength g for r+ = 1, 
r_ = and temperatures T = 5A (left) and T = 0.05A (right). The labelling of the curves is 
the same as in figure [2] For high temperatures (left), compared with 7 a/ 3, the temperature T 
simply rescales the D a p-ax\s. In the low temperature regime (right), the curves corresponding 
to sub-ohmic dissipation can present a minimum, while for ohmic dissipation the initially 
monotonically increasing curve is turned into a monotonically decreasing one. In the case of 
super-ohmic dissipation, the influence of temperature is smaller than for sub-ohmic and ohmic 
dissipation. 



is depicted). For low temperatures there exist corrections to (fTTj) : In the case of 
sub-ohmic dissipation, D aa (D a p, a 7^ 0) still diverges to +00 (—00), but now 
a local minimum (maximum) can occur. For ohmic dissipation, decreasing tem- 
perature turns the initially stricly monotonically increasing (decreasing) curve 
of D aa (D a p) into a strictly monotonically decreasing (increasing) one. For 
super-ohmic dissipation the situation is more complex, but in general the tem- 
peratures necessary to cause deviations from the high temperature behaviour 
are lower than those for subohmic and ohmic dissipation. Independently of the 
temperature, D aa is always positive, while this need not to be the case for D a p 
(a ^ j8). As one sees, the coupling of the oscillators leads to a rich behaviour 
already at the level of the coefficients jap and D a p. This richness arises from 
the effect of the high frequency modes of the environment on the system; on 
the contrary, 7 a/ g and D a @ are independent of the interaction strength g in the 
Caldeira-Leggett limit (which corresponds to a long wavelength approximation, 
see 437}). 

Let us now discuss the physical implications of the above mathematical 
observations. Because of the formal similarity of "f a /3 and D a p we will limit our- 
selves to the dissipation coefficients 7 Q( g . From a quantum-optical point of view, 
the fact that 7 Q( g and D a p are double Fourier transforms allows us to interpret 
them as a resonance of the environmental eigenmodes at the eigenfrequencies 

W r+ ± \J ' k 2 + r 2 " of the system. As the eigenfrequency r + — k 2 + r 2 _ is 

going to zero for g g^ , it is supposed that the divergency for g — ¥ g^ in the 
sub-ohmic case is connected with the divergence J'(uj —►()) = +00. If one keeps 
in mind that "f aa discribes the dissipation to the reservoir of the oscillator a, it 
is clear that j aa is always positive, for energy is transferred to the heat bath 
until the oscillator is slowed down to zero. That 7 Q( 3 (a / which describes 
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Figure 4: Dissipation in the case of distinct reservoirs. The oscillators dissipate energy to 
their heat bath at the rates 711 and 722, respectively, while due to the coupling, energy is 
exchanged between the oscillators at equal rates 712 = 721 • 



the dissipation of oscillator a due to the coupling to oscillator /3, can become 
non-positive is not a problem, since in the two-oscillator system 7 a( g = 7^ a and 
an equal amount of energy is exchanged between the oscillators. This situation 
is illustrated in figure |U 

When both oscillators are coupled to the same heat bath, the situation is 
different: As 712 = 721 can be negative, it is possible that "fi = 7x1 + 721 and 
72 = 722 + 712 can also become negative. Adding (joT)) and one finds out 
that the possibility of 71 and 72 to become negative depends on the sign of r_ . 
More precisely, the only term in the expression for 71 that can be negative is 

W 1 k 2 + r^ — K—r^, which is only the case for r_ > (corresponding to Qi > fig). 

In the same manner, the only negative term in 72 is yj k 2 + rl — K+r_ for r_ < 
(f2i < SI2)- For r_ = (f2i = f^) the coefficients 71 and 72 are always positive. 
Hence, 71 and 72 cannot become both non-positive. Moreover, 71 + 72 is always 
positive, which ensures that 72 > I71I (for 71 < 0) and 71 > I72I (for 72 < 0) are 
hold, i. e. energy is lost to the reservoir at a higher rate than it is gained. Under 
special circumstances (that depend in a complex way on the parameters s, n 
and r±) 71 (or 72) is indeed negative (figure [SJ. Let us assume that fii > 
and 71 < 0. In this case, oscillator 1 gains energy from the common heat bath 
at the rate |7i| (figure At the same time, oscillator 2 looses energy to the 
heat bath at a bigger rate 72 > |7i|. Due to the coupling of the oscillators, an 
excess of energy is transferred from oscillator 1 to oscillator 2. This unexpected 
mechanism leads to the effective dissipation of both oscillators and is only pos- 
sible because the oscillators are coupled to the same reservoir. Evidently, such 
a situation cannot be found by using the Caldeira-Leggett approximation, for 
which the coefficients "/ap are always positive or zero. The same arguments hold 
in the case of the decoherence coefficients D a/ 3. 
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Figure 5: Dissipation coefficients 71 (left) and 72 (right) for a common reservoir (here: r_|_ = 5, 
r_ = 3). The curves are plotted for different dissipation exponents s: dotted s = 0.5, 
continuous s = 1, dashed s = 1.25, dotted-dashed s = 1.5. For some values of g the coefficient 
71 can be negative, while 72 is always positive for r_ > 0. 
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Figure 6: Dissipation in the case of a common reservoir. For r_ > and special values of r±, 
s and k it is 71 < 0, i. e. oscillator 1 gains energy from the common reservoir. Nevertheless, 
71+72 > ensures that the other oscillator looses energy to the heat bath at a higher rate 
72 > |7i|- The coupling between the oscillators leads to the transfer of energy from oscillator 
1 to oscillator 2, so that on the whole both oscillators are effectively subject to dissipation. 



4-2. Reproduction of the Classical Motion and Quantum Corrections 

In this subsection, we will solve the differential equation system (|38[) for the 
expectation values of position for ohmic dissipation. This can be done ana- 
lytically by using standard methods, but the resulting expressions are lengthy. 
We have therefore refrained from stating them explicitly. Instead of the above 
exponential cut-off we will use a Lorentz cut-off function as in (|36[) . which has 
the advantage that there exist closed expressions also for fij^g and / a/ g (only for 
the case T = A/(2tt)) 1 and not only for 7 Q ^ and D a p . Note that the results of 
subsection 4.1 for s = 1 can also be obtained by using the spectral density (131))) . 
The choice of the cut-off function is therefore physically irrelevant. If one plots 
the solutions of (|3"5)) , one obtains the classically expected result: two coupled 
damped oscillations (figure [7]). In this example, due to the coupling, the oscil- 
lator with the lower frequency will follow the motion of the oscillator with the 
higher frequency. Figure [7] shows that the evolution of the mean position is not 
the same for distinct and common reservoirs, e. g. the low frequency oscillation 
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Figure 7: Temporal evolution of the expectation values of position for ohmic dissipation 
(A = 1, M = 1, Qi = 1, n 2 = 2). The initial conditions are (Xi)t=o = (^"2)4=0 = 1, 
{Pl)t=o = {P2)t—o = 0. The motion of oscillator 1 is depicted in black, the one of oscillator 2 
in red; continuous curves represent distinct, dashed ones common reservoirs. On the left: two 
coupled oscillators (g = 0.5, 70 = 0.2) which show the classically expected damped oscillations. 
The difference between distinct and common reservoirs is most striking for vanishing coupling 
g = (on the right, 70 = 0.5): For distinct reservoirs one obtains two independent damped 
oscillations, whereas in the case of a common reservoir the two oscillations are coupled by a 
reservoir-induced interaction. 



takes longer to follow the high frequency oscillation in the case of a common 
reservoir. The differences between the two heat bath models are the bigger, the 
greater the parameters 70, A and g are. As already mentioned in section [3j the 
case of vanishing coupling g — is particularly interesting: For distinct reser- 
voirs, one obtains two independent damped oscillations, whereas in the case of 
a common reservoir a reservoir-induced interaction couples the two oscillations 
(figure [7]). Before concluding this section, it has to be said that (pIS)) alone, with 
arbitrary coefficients f2^o and 7 Q ^, also admits a bunch of unphysical solutions 
(e. g. exponentially increasing oscillations). Only a physically reasonable spec- 
tral density, which enters in the quantum mechanical calculation of f2^o and 
jap , ensures physical solutions (for the calculation of the spectral density from 



microscopic models see |23[). 

To conclude this section, we briefly discuss how the interaction between the 
oscillators influences the quantum mechanical properties of the system. For 
this purpose, we calculate the quantum mechanical uncertainties AX a (t) = 
{X^)t — (Aa)t from the reduced density matrix for two initial Gaussian wave 
packets, described by the density matrix 

fe(x,x / ,0)= TT , 1 expf- X } +X " \ (72) 

As in the case of one harmonic oscillator [22j, the uncertainties saturate after 
some time. But unlike the expectation values (X) t , which end in the same fi- 
nal state of vanishing amplitude, no matter how strong the interaction between 
the oscillators is, the uncertainties reach different values for t — > 00, depending 
on the interaction strength (figure [3 these values do not depend on the initial 
conditions) . This means that quantum mechanically an effect of the interaction 
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Figure 8: Temporal evolution of the uncertainties of two coupled oscillators in separated 
reservoirs for ohmic dissipation (7 = 0.2, A = 1, T = l/(27r), M = 1, fii = 1, f22 = 
2). The initial uncertainties are chosen to be AXi(0) = AX2(0) = 1. Continuous curves 
represent AXi, dashed ones AX2. The saturation values are not the same for uncoupled 
(black) and coupled (red, g = 1) oscillators. For small times f the uncertainties calculated by 
using the Caldeira-Leggett approximation (green, g = 1, see J 371) ^1 differ not much from those 
calculated with the help of the Born-Markov approximation. For larger times the differences 
are considerable. 

remains even at t — > 00. In figure [5J one also sees that for low temperatures 
the uncertainties obtained by the Caldeira-Leggett approximation differ consid- 
erably from those calculated by using the Born-Markov approximation, espe- 
cially for large times. It is also worth noting that for the studied temperature 
T = A/(2-7r) the Caldeira-Leggett uncertainties systematically underestimate 
the Born-Markov uncertainties and relax more rapidly into their final values. 

5. Summary and Outlook 

In the present paper we worked out how to calculate the density matrix 
for a system of linearly coupled oscillators in a bosonic heat bath using the 
Born-Markov approximation. By having done this, we complemented previous 
results from [19J, which are only valid in the Caldeira-Leggett limit. As long 
as the Born-Markov approximation holds, we are now able to access the non- 
ohmic, low temperature regime. In the case of the two-oscillator system we 
found that the expectation values of position reproduce the classically expected 
result of coupled damped oscillations with certain renormalized coefficients. As 
for the quantum uncertainties, they are modified by the coupling even in the 
limit t — y 00. Moreover, we have explictly seen that the coupling to a common 
reservoir leads to reservoir-induced interactions. In this context, using a specific 
spectral density, we found an interesting dissipation mechanism for a common 
reservoir, which is based on the fact that energy can be tranferred via the heat 
bath. This mechanism is due to the influence of the high frequency modes 
of the environment on the coupled system and is therefore not observable in 
the Caldeira-Leggett limit. The above results should encourage further studies 
in decoherence of interacting systems, in particular regarding the coupling to 
a common reservoir 0, We note that the presence of several reservoirs 
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leads to interesting effects already for classical particles [24J . The calculation of 
other interesting physical quantities, like the entanglement of the subsystems, 
is postponed to a later work. In addition, it would also be interesting to study 
the effects that occur in the strong coupling, non-Markovian regime. 
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